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I. Introduction. 

Instabilities of extended systems lead to a large number of interesting phenomena and 
in particular to the appearance of well defined structures^ 5 ' 12 ). We will consider in this 
paper the case of discrete time iteration of lattice systems. These are dynamical systems 
defined on the phase /°°(Z d , R) which will be abbreviated below by The time evolution 
is given by a map $ of l°° into itself for which we will make the following basic hypothesis 
which also includes the dependence on a real parameter. 

HI $ : Rxl°° — ► l°° . $ is C 2 in its first argument and C 4 in the second one on some 
neighborhood of the origin. Moreover $ commutes with the translations on Z d . 

The time evolution for a fixed parameter rj is given as follows. If u 1 G l°° (t G N) is 
the state of the system at time t, the state of the system at time t + 1 is given by 

yt+i _ $(77^*) . 

There are many examples of such dynamical systems. We will just mention the class 
of coupled lattice maps and also notice that many algorithms used for solving numerically 
partial differential equations are of this form (for example finite difference schemes). 

One would like of course to understand the properties of this dynamical system and 
in particular the large time dynamics. In finite dimensional dynamical systems, varying a 
parameter is a very efficient way to discuss the increase of complexity of the time asymptotic 
dynamics. This increase of complexity can occur at well defined threshold values of the 
parameter called bifurcations for which detailed results are available in particular for the 
bifurcation of stationary solutions^ 8 ) . An important tool in bifurcation theory of stationary 
solutions is the normal form equation which provides a model for the dynamics. Moreover 
for generic systems these normal forms are in some sense universal. 

For spatially extended continuous systems which depend on a parameter, amplitude 
equations have been derived in various situations^ 5 ' 12 ), and rigorous derivations of the 
Ginzburg-Landau equation in one space dimension have been obtained*- 3 ' 15 ' 11 ' 13 ' 14 ). The 
result is similar to what occurs for normal forms in finite dimensional dynamical systems, 
namely the space-time function reconstructed using the solution of the amplitude equation 
stays close to the solution of the true evolution (with the same initial condition) for a large 
time. However this method usually provides only an approximation of the true solutions 
and for very large times, neglected small errors usually grow due to instabilities^ 10 ). Never- 
theless, just before the approximation becomes too inaccurate one can restart the process 
with a new initial condition which is the present true solution. This provides a shadowing 
of the true solution by a sequence of chunks of solutions of the amplitude equation. We 
refer to the previously mentioned references for more detailed discussions of these results. 

The shadowing of a solution of the true equation by a reconstructed function from a 
solution of the amplitude equation is only valid for initial conditions of a rather limited 
form, namely those which can be reconstructed from an initial condition of the amplitude 
equation. Eckhaus'- 6 ) has shown that this condition can be somewhat relaxed and that 
after a relatively short time any small enough initial condition leads to a solution having 
the right form. 

We will consider below all the previously mentioned results for the case of discrete 
dynamics of one dimensional lattice systems. We will recover all the results obtained in 
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the continuous case and extend them in various directions. In order to do this we will use 
a renormalization group approach^ 1 ' 7 ). Besides its systematic simplicity this method will 
also have the advantage of providing a rigorous basis for the universality of the amplitude 
equations. In the finite dimensional situation transversality arguments can be used, how- 
ever in the infinite dimensional case the equivalent results are not known, and universality 
provided by the renormalization group method is a good replacement. 

Another advantage of the renormalization group argument is that it will allow us 
to deal with initial conditions which are small but of order unity, independently of the 
smallness of the parameter. As we will see below the renormalization group flow draws 
an initial data of size of order one to a configuration of size ^/fj. At the same time the 
renormalization group flow puts the solution in a resonant form providing a natural proof 
of the previously mentioned Eckhaus result. Finally when the renormalization group stops 
improving the size and form of the solution we get the amplitude equation together with 
a rigorous proof of the associated shadowing properties. The same ideas can probably be 
applied in the continuous space and time case and will eventually provide some improve- 
ments in term of simplicity and naturalness of proof as well as better bounds. Different 
approaches to amplitude equations using renormalization have been discussed by Golden- 
feld et al/ 4 -*. However since there is already a large literature for the continuous case we 
will not consider this situation. Note also that we will treat a rather large class of systems 
whose continuous equivalent would include pseudo differential operators. 

We will mostly deal below with one dimensional systems (d = 1) although the same 
method works in higher dimensions, at least when the amplitude equation is not too 
involved. Systems of coupled equations (corresponding to vector valued fields) can be tee 
treated similarly. 

Since we will be interested in instabilities of stationary solutions, we will now assume 
for simplicity that the zero element of Z°° is a stationary solution of the dynamical system. 
In many cases this can be realized by a translation in function space. 
H2 $(77, 0) = for any 17 e R. 

We will make some assumptions on the linearized evolution around this fixed point. 
H3 We will assume that the operator L = D 2 &(0, 0) which is bounded in l°° and com- 
mutes with translations can be represented as the convolution with a (real) sequence I El 1 . 
Let 1(9) be the Fourier transform of I 



1(6) = l ne 

n6Z 



ind 



We will also assume that 



nez 



n| 10 |/ n | < 00 , 



and that \l(0)\ < 1 with \l(u>)\ = 1 for a unique uj e]0,7r[. We will also assume that 
l(uj) 7^ —1 and uj £ {tv/2 : 3tt/3}. Moreover, we will make the non degeneracy assumption 



D = 



1"(oj) _ l'(uf 
l(u) l(u) 2 



7^0 
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As we will see later on, this last number is always non-negative. Notice that 1(oj) = 1(—oj). 
The particular cases which are excluded in the above hypothesis correspond to some strong 
resonances. We will not consider them further below although some interesting results can 
probably be derived in these situations using the ideas of the present paper. 

We now need some hypothesis on the higher order derivatives of 
H4 The translation invariant bilinear and trilinear operators Q = D^iO, 0)/2 and C = 
D|$(0, 0)/6 are given by a double (respectively triple) real sequence (q n ,m) / 1 (Z 2 ) 
(respectively (c n>m)P ) e l 1 ^ 3 )), namely 

Q(u)r ^ ^ Qr—n,i — m^n^m and C(tt) r ^ ^ C, — Ujl — m ^ — pU n U m Up . 

(n,m)eZ (n,m,p)€Z 

We will assume that these sequences satisfy 

£ \q n ,m\(l + H + |m|) 6 + £ \c n ,m, P \(l + \n\ + \m\ + Ip]) 4 " < oo . 

(n,m)£Z (n,m,p)gZ 

We will also make the non linear stability assumption 

= _ / c(a^, -cj) 4g(o;,0)g(a;, -oj , 2g(-o;, 2o;)g(a;, u) \ 
9 V *M 3Z(c^)(|Z(c^)| 2 - Z(0)) + 3Z(w)(Z(u;) 2 -Z(2u;)),/ ' 

where 

c(0i, 02, 3 ) = £ c p , g , r e- l W 1+ ^+^) and g(0x, 2 ) = £ g^e^i-H*.) . 

Finally we will also require that the dependence on the parameter allows the instability 
to occur. 

H5. We will assume that the operator L\ = DiZ^^O, 0) which is bounded in l°° and 
commutes with translations can be represented as the convolution with a (real) sequence 
m G I 1 . We will assume that 

|n| 3 |m n | < oo . 

nez 

Let m{9) be the Fourier transform of m, we will also assume that 

1{0J) 



Note that as for the function / we have m(—9) = m(6). 

A simple example satisfying the above hypothesis can be satisfied as follows. Let / be 
a bi-infinite sequence of real numbers satisfying H3, we define the map $ by 



®(V,M)n = VUn + ^jln-qUq - (u n f . (1 .1) 
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In this particular situation, HI and H2 are easily verified. 

This is a natural generalization of the discrete Swift-Hohenberg evolution which is 
given by 

< +1 = (1 + VX - P((l + 7A) V) n - «) 3 (DSH) 
where A is the discrete Laplacian 

(Am)„ = u n+1 + u n -x - 2u n . 

In order to explain some of the arguments in a simpler situation we will also sometimes 
make the following hypothesis. 
H3' / is even and satisfies H3. 

Note that in this case 1(9) is even and real, 1(±uj) = 1 and V(±uj) = 0. 
For the particular case of (DSH) one has 

l((f)) = 1-^(1-27 + 27 cos (p) 2 , 

and the hypothesis H1-H5 and H3' are satisfied if one assumes 

7 > 1/4 and 2/p > max(l, (1 - 4~f) 2 ) . 

We now formulate the main result of this paper for the evolution (1.1). 

Theorem 1.1. Assume H3', then there are constants R > 0, 1 > T]o > 0, c\ > 0, • ■ • , c 4 > 
such that for any rj e]0, rjo[, for any u E l°° such that \\u\\i°° < R, if (M*)teN denotes the 
orbit of u under the discrete time evolution (1.1), then 

i) there is a positive number T = T(u, rj) such that for any integer t > T 

11^11;°° < cin 1/2 . 

ii) For any integer t > T there is a function A(r, x) from R + x R to C satisfying 

d T A = d 2 A + A — A\A\ 2 (GL) 
such that for any integer s with < s < C2?? _1 log?] -1 we have 

\\u t+8 -v s \\i°o < c 3 ri 1/2+Ci . 

where 

v s n = 3" 1 / 2 e^ n ?7 1/2 A(?7s, ^^n^l/D) + c.c. 

Theorem (LI) can be understood in the light of the basic principles of bifurcation 
theory. First of all we have an "unstable" direction in phase space which is the sequence 
(e lujn ) (and its complex conjugate). As in finite dimensional bifurcation theory, the non 
trivial solutions are sought in the form u n = A n e lujn , where A is a slowly varying function 
of space and time (it is a non trivial consequence of our theorem that locally there are 
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no other long time solutions). The equation for A should ensure that in the evolution 
equation for u all resonant terms (i.e. terms containing e ±MiJn ) disappear. These resonant 
terms come of course from the linear part of the evolution equation but also from the non 
linear part. If we examine more carefully the second contribution, we see that quadratic 
terms in u (and all even higher order terms) should not contribute, whereas cubic terms 
(and higher odd powers) generate resonant terms. If the solution is small, cubic terms 
dominate. If we assume that the dispersion relation is non degenerate at the critical wave 
number (here D > 0), one is naturally lead by scaling arguments to the equation 

d t A = A" - A\A\ 2 . (F.P.) 

Note that this equation has the one parameter family of scale invariance x — > \x, t — > A 2 t, 
A — > XA. This equation can be considered as a fixed point of a renormalization operation. 
An important property of (F.P.) is that all solutions with bounded initial conditions tend 
to zero when t tends to infinity. A term a A added to (F.P.) corresponds to a relevant 
direction^ 1 -'. This is how the complete Ginzburg-Landau equation appears. The proof of 
Theorem (LI) is essentially a control of the R.G. flow in the vicinity of (F.P.) done from 
the point of view of the solutions i.e. proving that the unstable manifold is transversally 
stable. 

We implement these ideas in section II. We first give estimates on various operators 
which are needed to control the time evolution of the remainder. The recursive renormal- 
ization argument is given at the end of the section. 

In section III we will give the equivalent result under the more general hypothesis Hl- 
H5. Technically the proof is not much different (except for the treatment of the quadratic 
term). However one has to deal simultaneously with several details which somewhat ob- 
scure the main ideas, this is why we have chosen to explain first in section II the simple 
case, and explain in section III the necessary modifications for the general case. 
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II. Renormalization proof of Theorem II. 1. 

In this section we will prove Theorem 1.1. We will therefore always assume that 77 is 
a positive number. 

We start by fixing a renormalization factor S which is a positive number larger than 
one. Although the final result will be essentially independent of S, it is convenient in the 
present discrete space-time setting to take for S and integer (for example S = 2). This 
avoids using integer parts of various time scales. We can now define a sequence of scales 
(5 n ) nfE N for the wave numbers by 

In other words we have a sequence of space scales given by (5~ 1 ) n6 N. As we will see below, 
since we have a "trivial" scaling, this induces a a sequence of time scales (5~ 2 ) n6 N, and 
a sequence of field scales given by (5 n ) n6 N • We will also need below a sequence of scales 
(#n)neN for a window of unstable wave numbers. It is convenient to choose this sequence 
of the form 

where <; is a positive number smaller than one and small enough. We will for definiteness 
choose q = 1/5. 

As usual with renormalization group arguments, all the steps are similar except for a 
change of scales. Therefore we will start by assuming that a scale has been fixed (i.e. an 
n in the previous sequences). Hence we have a scale 5 = 5 n for the solution, namely we 
will assume from now on that the initial condition u satisfies 

S~ x 5 < < ^ • 

We also have a scale 9 for the normalized width of the window of wave numbers. One 
should think of 5 as a small number (but possibly of order one, i.e. a negative power of 
and 9 as a large number such that 59 4 is small (i.e. 9 = 5 1 / 5 ). Explicit constraints 
on 5 and 9 will be given later on. In particular we will always assume that 59 < 1/8 (see 
appendix B). 

We now introduce three operators whose sum is the identity and which are essentially 
projections on sequences with wave numbers near u. They are constructed as follows. 
First we consider a cut-off function if; which is infinitely differentiable, non negative with 
support in [—2,2] and which is equal to one on the interval [—1,1]. We then define a 
sequence Kg^g by 

K s ,e(n) = 59FiJj(n59) , 

where T denotes Fourier transform on the real line. Convolution by Ks,e localizes the wave 
numbers in an interval of length 469 centered at the origin. In particular, for 59 < 7r/2, 
since that tp has compact support we have for — tt < 4> < tv 

Y,K s , e {n)e in * = *1>{9- 1 5- 1 <I>) . 

n 

We now define three operators Pg~g, Pg~ e and R$ : g by their kernels, namely 
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P s ~ e (n, m) = e-^ n K s , e (n - m)e^ m , 

and 

Rs,e = I - P£ e - Pf j0 ■ 

From the above remark about the action of the convolution by Ks } g, we see that P$ e 
localizes the wave numbers around u> and Pf around — u. In particular, if 5 is small 
enough, we obtain two disjoint intervals of wave numbers. 
We now define a constant C\ by 



Ci = 2 1+ sup ||a^vili • 

V 0<j<5 / 

This constant will be a convenient bound for several quantities below. 

Lemma II. 1. For any < 56 < 7r/2, the operators Pg~ e , Pg~ e and Rs,e have a norm 
bounded by C\ in 

The proof is an immediate consequence of the fact that the function if; belongs to the 
Schwarz space S. 

We now decompose a sequence u G l°° using the above three operators. 
Lemma IL2. For any real sequence uEl°°, there is a function B e C 4 (R, C) such that 

e-^ n P+ 9 (u) n = B(5n) \/n G Z , 

and satisfying 

max{0- J '||d*£|| I oo ;j = 0,---,4} < Ci||u||ioo . 

We have also 

e lLJn P^ e (u)n = B{5n) \/n G Z . 

This is an immediate consequence of the definition of P^ B and interpolation. In order 
to fix entirely B one can use the following explicit choice 



B(x) =56^2^ ( xd ~ mSd ) e ~ 



■imu> 
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We also define the remainder r° by 

r° = Rs,9U . 

We now denote by A the solution of the Ginzburg Landau equation 

d T A = A" + t]5~ 2 A — A\A\ 2 , (II A) 

with initial data 

A(0, x) = V3 S^Bixy^D/^) . (II A') 
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It will be convenient in the sequel to denote by J± the sequence 

A* n = 3~ 1/2 e inuj A(5 2 t, n5^2/D) + c.c. (II. .2) 
Finally, the equation for the evolution of the remainder term 

r* = u* - 5 A 1 

is given by 

r' n +1 = (Lr% - 35 2 r* K) 2 - 3(r* ) 2 A\ - Kf + ft* , (JJ.3) 
where we have denoted by I the operator 

Iv = I *v+ rjv , 

and the sequence TZ 1 is given by 

K = S(LA% - d 3 (A t n ) 3 - SA^ 1 . (II A) 

Using equation (II. 1), we have immediately the following estimate for the forcing term 

Lemma II. 3. There is a constant C 2 > C\ such that for any 56 4 < 1, and 5 > n 1 / 2 , if B 
is a function of x satisfying 

max{e- J \\d J B\\ la o ■j = 0,---,A}<6C 1 , 

if A denotes the solution of (II. 1) with initial condition (III'), and if A 1 is constructed 
from A as in (II. 2), then for any t G N we have 

11^*111°° < C 2 , < C 2 5 3 , ||S* + C*||i=c < c 2 5 4 e* , 

where the sequence C* is given by 

C f n = 3~ 3/2 5 3 e 3inLJ A 3 (5 2 t, n5^2/D) + c.c. 

Proof. It follows easily from (II. 1') that 

\\A(0, -)||/- < VSC! . 

The Lemma follows then at once from Lemmas A. 3 and A. 4 of appendix A using 56 4 < 1. 
We now start estimating various operators needed below in the proof of Theorem 1.1. 

Lemma II. 4. There is a constant C3 > e such that for any < 77 < 1 and for any t G N 

we have 

\\L%oo < C 3 e 7 ' (t - ,? " 1) . 



The proof is given in appendix B. 
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Lemma II.5. For any t eN, < 77 < 1 and 1 > 5 2 9 2 > 4Dn we have 

WL'RsMl <C 4 e~ tDs2e2 / 8 . 



The proof is given in appendix B. 

It will be convenient to denote by L t the linear operator given by 



(L t v) n = -SS'iAiYvn - SSrlA^Vn - {rlYvn . 

Note that 

r* +1 = (LrX + (L t r% + TZ l n . (I/.5) 
Lemma II. 6. We have for any t > 

\\L t \\i°° <5(5 2 ||^ t ||f 00 + ||r*||jL) • 

This is obvious from the definition of L t . 
Lemma II. 7. If < t < t < rj~ l we have 



\\(L + L t )-.-(L + L t0 )\\ l -<C 3 (l+ sup 5C 3 (<5 2 || 2 4 s ||f 00 + ||ri^)') 

V t <s<t J 



Proof. We first expand the product (L + L t ) • ■ ■ (L + L to ). It is a sum of 2* to+1 terms of 
the form 

where the numbers (. are equal to zero or one. Using Lemmas II. 4 and II. 6 it easily follows 
that 

\\L 1 -<*Li t ...L 1 -<*L%\\ l -<Ct 1 5 k ( sup (5 2 ||^ s || 2 00 + ||r s ||^)') 

\t <s<t J 

where k is the number of (. equal to one. Therefore 

|| (L + L t ) ■••(!/ + LJ\\^ < |5 ( t ~ t ° k + 1 ) C k 3 +1 5 k ^supJ6 2 \\A s \\^ + llrl?-))* 

and the result follows. 

lit — to > we can split this time interval in various sub-intervals of length [i] -1 ] 
plus a remainder. We immediately obtain the following corollary. 
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Corollary II. 8. For any < to < t and rj < 1/2 we have 

\\(L + L t )-..(L + L t0 )\\t~ < <:I7 i+2^( t -*o) e 503(*-*o+i) ^P.o^^.C^-ll^llf^+llz:-!!?^) 



We now start estimating the remainder r*. This will be done in two steps. Namely, 
during a time of order 5~ 2 we do not loose too much on the estimate of the remainder. 
After this time has elapsed and up to a time of order 0(1)5~ 2 log 5 _1 the estimate improves 
because the remainder is in a contracting subspace of the linear evolution. When we reach 
this time (after which we may somewhat loose again on the estimate), we perform a 
renormalization step and start the whole process again. 

Lemma II. 9. There is a constant 1 > C 4 > such that if C 4 > 5 2 > 77, if t < C^5~ 2 
and 6 < 5" 1 /4/(2C 2 ) then 

||r*||/oo < 4C 2 C 3 d . 



Proof. Let the non decreasing sequence of positive numbers a t be defined by 

a t = 5~ x sup Hz^llz 00 . 

0<s<t 

Note that by definition of C\ we have cr < C\ < 4C 2 C 3 . We have from equation (II. 5) 
the identity 

t 

r t+1 = (L + L t ) ■ ■ ■ (L + L )r° + ^(L + L t ) ■ ■ ■ (L + L j+1 )K j . 

3=0 

Using Lemma II. 3 and Lemma II. 7 we obtain 



\r t+1 \\l< 



t 

< C 3 (1 + 5C 3 5 2 (C 2 + a 2 )) t+1 ||r || ;oo + ^C 3 (l + 5C 3 5 2 (C 2 + a 2 ))'' 3 ||S J "IU«= 

3=0 

since t5 2 < 1. The result now follows recursively from the choice of the constant 

log 2 



C 4 



170CfC 2 2 ' 



We now improve the previous result using the support properties of the spectrum of 
r° and of TZ 1 . We first derive an equivalent of Corollary II. 7. 
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Lemma 11.10. There is a constant C5 > 1 such that if r\ < 1/2, 59 < 1, 5 2 > rj, and t is 
such that for < s < t we have \\r s \\ < AC2C35, then 

|| (L + L t ) • • • (L + L )Rs,o\\io° < C 5 9- 2 e c ' 5ts2 + C^- uH ^ c * . 



Proof. We first observe that 

t 

(L + L t )---(L + L ) = L t+1 + J2(L + Lt) ■ ■ ■ (L + L S+1 )L S L S , 

s=0 

where as usual, the elements with indices out of range are equal to identity. We now apply 
Lemmas II.5-II.7 and Corollary II.8 to get 

|| (L + L t ) ■ ■ ■ (L + L o )R 5y0 IIjoc < C±e- tm2B2 / 8 + 5(C 2 2 + 16C 2 2 Cf )5 2 
t-i 

^ e -sD5 2 9 2 /8 C 1+Mts) e 5C 3 (t-s)5 2 (CZ + 16C%Cl) 
s=0 

from which the result follows immediately. 

Corollary 11.11. There are constants Cq > C5 and C? > such that if 5 is small 
enough so that the numbers ±3uj (mod 2ir) do not belong to [u> — 265, u> + 295] U [— u> — 
295,— uj + 295], then if t is such that for < s < t we have \\r s \\ < 4C2C35, and if 
C 7 5~ 2 \og5~ 1 > t > C A 5~ 2 we have 

WrH < U59- 1 ' 2 . 



The proof is similar to the proof of Lemma II. 9 using Lemma 11.10, Corollary CI, and 
the relation 

t 

J2(L + L t )-..(L + L j+1 )TZ j = 
3=0 

t t t 

Y.L^W+Y, E (L + L t )---(L + L a+1 )L a L-*- 1 W . 

j=Q j=0 s=j+l 

We can now give the recursive proof of Theorem 1.1. We start by selecting an integer 
no large enough so that all the constraints on 5 = 5 n and 9 = 9 n in the previous Lemmas 
are satisfied for n > no. We also assume no large enough so that log^g -1 > (S + l) 2 + 1. 
We define R = 5 . Note that the choice of n is independent of n. We now define r]o by 
Vo = S-%/8. 

Since the proof is recursive we will therefore assume that we have reached a situation 
corresponding to scale 5 n (n < no). In other words we assume that u° satisfies 

\\u < 5 n . 
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We will also assume that 

||u°||joo > S^dn , 

otherwise the choice of n is not optimal. Finally since we are for the moment interested in 
the first part of the Theorem, we will also assume that S 2 > r\. 

We now consider the time evolution starting from this initial data yr, and we will 
apply the previous Lemmas with 5 = 5 n and 9 = n = We first apply Lemma II. 2 to 
get a function B with the corresponding estimates on the derivatives, and we define the 
function A(t,x) solution of (II. 1) with initial data (II. 1'). We also define the remainder 
r°. 

The time evolution will now be split into two parts. We define a first time 

tn = [C4<5 n 2 ] . 

On the time interval [0,£ n ], we have from Lemma II. 9 

llrijoo <AC 2 C 3 S n . 

We now define a second time by 

T n = [C 8 (5- 2 log(5- 1 ] , 

where CV > C% > and C% is small enough so that for t < C%5~ 2 log <5 _1 , all the ex- 
ponentially growing factors which appear in all the previous Lemmas are smaller than 

On the time interval [t n ,T n ], we can now apply Lemma 11.11. In particular, there is 
a time T n which satisfies 

t < T < T 

where we have 

||r r "|| ; oo < S- 1 5 n /2 , 

and also using Lemma A. 2 

p-^ 2 A{8lT n , Oik- <S~ x /2. 

Therefore 

and we can now repeat the argument with n + 1. 

The proof of the second part of Theorem 1.1 is essentially similar. We first define an 
integer iV as the smallest integer such that Sn < Srj 1 / 2 . We can now repeat the above 
proof and the result follows since Tn is of the right order. 
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III. The General Case. 

In this section we will prove the following theorem which is a generalization of Theorem 
1.1. 

Theorem III.l. Assume hypothesis H1-H5. Then there exists constants R > 0, i]q > 0, 
Ci > 0, • • • , C4 > such that for any rj e]0, t]o[, for any u G l°° such that < R, if 

{u t )te~N denotes the orbit of u under the discrete time evolution $, then 

i) there is a positive number T = T(u, rj) such that for any integer t > T 

< cir] 1/2 . 

ii) Define the constants a, f3, a, 7 and V by 

Dl(u) pl\^) ly^) 

and 

n_ ^ f c(uj,uj,-uj) 4q(u,0)q(uj, -u) , 2q(-uJ,2uj)q(uJ,uj) \ 
f " 9Up) 3^M(|/M| 2 -/(0)) + 3^(a;)(;M 2 -/(2 W ))y' • 

For any integer t > T there is a function A(t, x) from R + x R to C satisfying 

d T A = (1 + io^A + A - (1 + z/3)A|A| 2 (CGL) 

such that for any integer s with < s < C2?7 _1 log?] -1 we have 

where 

V s n = jfrt+rnt+utn) ^1/2 /Z^^, 7] 1 / 2 y^/D(n + Vt)) + CC . 

We recall that the constants D and p have been defined in hypothesis H3 and H5 
respectively. 

Lemma III. 2. The number D is non negative, and the number V is real. 

Proof. Let Q(9) = l(6)/l(u>), this function is equal to one in 9 = uj and this is a point 
where its modulus is maximum. Let (jt and £j denote the real and imaginary parts of £. 
Writing that the first derivative of the square of the modulus of Q in u> is equal to zero and 
the second derivative is non-positive, one gets immediately (using (i(uj) = 0) 

c^m = o, and c^m + (c;m) 2 <o, 

and the lemma follows. 

From the differentiability of $ we can write 

$(77, u) = L'u + r]L"u + Q(u, u) + C(u, u, u) + G(r], u) , (^-1) 

where G is a nonlinear map, differentiable which is at least of degree 4 in u, at least of 
degree 2 in 77 and at least of degree one in rj(u) 2 . 

Our first goal is to eliminate the quadratic term Q. As we will see, this is not entirely 
possible, but it will be enough to eliminate the dominant part of this map. The main idea is 
of course that Q should not give rise to a resonant term, and therefore it should be possible 
to eliminate the quadratic terms by a nonlinear change of variables. An implementation 
of this idea is given by the following result. 
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Theorem III. 3. There is a (differentiable) quadratic operator S in l°° such that if S 
denotes the operator 

S (u) = u + S{u) , 

then on some neighborhood of zero S is invertible (with differentiable inverse), and the 
maps 

$(r},-) =S -1 o%.)o5 

satisfy the hypothesis H1-H5 with a quadratic part Q = -D^^KO, 0)/2 such that its Fourier 
transform q{<pi,<p 2 ) is identically zero on a neighborhood of the points (±u, ±u) and of 
the lines <pi -\- <p 2 = Q satisfies also 

+ \n\ + \m\) 6 \ 

Qn,m\ < 00 • 

n,m 

Moreover, we have 

D 2 £(0,0) = L', DiD 2 ^(0,0) =L" , 

and 

Dl$(0,0)(v,v,v)/6 = C{v,v,v) + 2Q{v, S{v)) , 

and the Fourier coefficients of this form decay also polynomially fast. 

Proof. By continuity, it follows from H3 that we can find a number h < uj/4 such that 
if (f>i and <p 2 belong to the intervals [±u — 2h, ±uj + 2h] then I {(pi + <p 2 ) — l{<pi)l{<p2) 7^ 0. 
Moreover, we can choose h small enough such that if \(pi + (p2 ± uj\ < 2h then either 
|0i ± u\ > 2h or \<p 2 ± u\ > 2h and l{(pi + <p 2 ) - l{(pi)l{(p 2 ) ^ 0. Let 

MO) =M0 + u J )h~ 1 /2) u)h~ 1 /2) , 

where ip is the cut-off function defined at the beginning of section 2. We now define a 
double sequence {T n>m ) by 

Since the function ip defined in section 2 is regular with compact support, it follows that 
the sequence (T n ^ m ) belongs to / 1 (Z 2 ). Note also that this sequence is real. Moreover, we 
have 

^|r n , m |(l + |n| + |m|) fc < 00 

n,m 

for any integer k. 

We now define the double sequence s by 

s = T * q , 

and the quadratic map S by 

S{v) n = ^2 S 



n—p,n—q'VpVq 



P,Q 
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The proof follows by easy computations. 

We will now drop the tilde on the various transformations, i.e. we will assume that 
the function q(<fii,<p2) satisfies the conclusions of the previous Theorem. We can recover 
the general result by applying the above Theorem. 

The proof of Theorem III.l is now very analogous to the proof of Theorem 1.1. It goes 
along the same steps of the renormalizaztion group using estimates analogous to those of 
Lemmas II.3-II.il. We will only indicate the non-elementary changes in the proofs without 
keeping a detailed account of the various constants, i.e. the notation 0(1) will be used for 
any constant independent of 5, rj and u. 

We first decompose u° as in section 2, namely we first apply Lemma II. 2 and we obtain 
a function B. As in section II, we set 

u* = 5 A 1 + r* 

where now 

A\ = e l(rTt+ujn \3g)- 1/2 A(5 2 t, 5^2/D(n + Vt)) , (1 1 1. 2) 

with A a solution of the complex Ginzburg Landau equation 

d T A = (1 + ia)A" + (p + i-f)r]5- 2 A - (1 + i(3)A\A\ 2 , (CGL) 

with initial condition 

A(0,x) = ^5~ 1 B(x^Dj2) . (CGL') 

The remainder r* satisfies the equation 

r t+1 = (C t +L t )r t + n\ (III. 3) 

where 

= SL'A* + vt^'A* + S 2 Q(A t ,A t ) + 5 3 C(A\ A\ A f ) + G( v , 6A*) - SA t+1 , 
and the linear operators L t and C t are defined by 

C t v = L'v + V L"v + Q(v, 25 A* + r*) 

and 

L t v = 3C(v, 8A\ 5 A 1 + r l ) + C(v, r\r l ) + [ d^D 2 G( V , 5 A 1 + . 

Jo 

We first describe the estimations on the linear operators. Note that Lemma II. 7, 
Corollary II. 8 and Lemma 11.10 follow as before if we have Lemmas II.4-II.6. 

The following equivalent of Lemma II. 6 is immediate from the above formulas. 

Lemma III.4. If ||<L4*||z°° + ||r*||joo < 1, then 

\\L t \\i°° <C»(l)(5 2 ||^||foc + ||r*||?oo). 



We now come to the equivalent of Lemma II. 4 and Lemma II. 5. 
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Lemma III. 5. There are two constants c > 1 and c' > such that for any r\ < 5 2 9 2 E 
]0, 6/2] (b is the constant of Lemma C3) and for any t > t > such that 

sup (5||^4*||/oc + ||r*||joo) < c'S , 

t <s<t 

we have 

HA •••£*„ III- <ce c "(*- t0 ) 

and 

||£ t -..A ^|| Zoo < ce - 52e2 ^)/ c . 

Proof. We first observe that C t = L + C' t where 

Lv = L'v + nL"v and C\ = Q(v, 25 A f + r*) . 
With this notation we have 

t 

£>t ■ ■ • A = ^* t0 + 1 + ^ ^* S ^s-l^s-2 • • • C to • 

Using Lemma C3 and Lemma II. 4 we get 

|| Ct ■ • ■ C tQ ||/oo < 

t 

O(l) e «»(*-*o+i) + ^ ^(^e-^^ll^ll^ + llr 8 - 1 !!^)!!^--.^!!^ . 

s=t + l 

The proof of the first part of the Lemma follows now recursively. 
To get the second part, we use the following formula 

C t ■ ■ ■ C t0 = L'-^ 1 

t-to+l 

+ 2^ L 't-n 1 lj L -t—n 1 -n 2 Ll ^ L 't-n 1 n k -k ^ 

k=l ™lH \-n k + 1 =t-t + l-k 

ni>0, •■■ ,nfc^!>0 

The second part of the Lemma follows now easily using Lemma II. 5 and Lemma C3. 

As explained before, Lemmas II. 7, Corollary II. 8 and Lemma 11.10 follow as before. 

We now come to the estimates concerning the forcing term TZ 1 . From Lemma C2, the 
term 5 2 Q(A\A t ) is at least of order 5 4 2 . For the term 11* - 5 2 Q(A t 1 A t ) we have the same 
estimates as in Lemma II. 3 and therefore we get bounds analogous to those of Lemmas 
II. 9 and 11.11. 

The proof of Theorem III.l then follows exactly the same steps as the proof of Theorem 
1.1 and we will not repeat them here. Finally one has to apply the inverse of the map S 
to conclude the proof. 
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Appendix A. Some estimates on the solutions of CGL. 

In this appendix we will derive estimates on the solution of equation (CGL). General 
estimates were derived before^ 2 -*, however the present situation is slightly better in the 
sense that the initial conditions will never be larger than some predefined constant. We 
will also restrict ourselves to dimensions one and two where the complex Ginzburg-Landau 
equation is known to have well behaved solutions for any values of the parameters a and 
(3. We will in fact obtain a stronger result than what we need in sections II and III. 
We consider the equation 

d t A = (1 + ia)AA + a(p + z 7 )A - (1 + i(3)A\A\ 2 , (Al) 

where a is a positive parameter smaller than a given positive constant S, and a, /?, 7 and 
p are fixed constants. Note that we could have as well rescaled a and fixed for example 
p = 1. We will not do this in order to have an immediate application of the result to 
section III. 

Lemma A.l. Let A{t, x) be a solution of (A.l). Assume that there is a number k > 2, a 
constant K > 1 and a constant 9 > 1 such that 

sup ^'11^(0, -)\\<K. 

0<j<k 

Then there is a constant M which depends only on K, S, a, f3, 7, p, such that for any 
t > 

sup 6- j \\d J 2 A(t, -)Hl~ < M . 

0<j<k 

Moreover, for any integer p < k/2 we have 

sup 9- 2j \\d{A(t, ■ < M . 

0<j<p 

There is a number U > which is independent of 9 and A such that for t > U the 
same estimates hold without the term 9~K We can take U as small as needed eventually 
increasing M. 

Proof. Let G t denote the function 

GJx) = , \, /0 e~ x2 / u . 

We then consider the integral equation 

A(t,-) = e iat G t *A(0,-) 

+ f dsG t . s * (a(p + n)A(s, •)-(! + 0)\A(s, -)\ 2 A(s, ■ )) . (A2) 
Jo 
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It is easy to verify that there is a number 1 > T > which depends only on k and K, S, a, 
[3, 7, p, such that this equation can be solved by the contraction mapping method on the 
time interval [0, T] with a solution in the ball of radius 2K in L°°. Moreover, the solution 
is also a solution of (A.l) with initial data A(0, ■ ). Taking the successive derivatives in x 
of (A.2) we get 

&A(t, • ) = e iat G t * &A(0, ■)+ f* G t . s * (a(p + ii)&A(s, • ) 

Jo 

-2(1 + ip)\A(s, ■)\ 2 d j A(s t -)-(l + ip)A(s, ■fd j A{s,-))ds + F j , 

where Fj is a combination of derivatives of A of lower order. Eventually taking a smaller 
T (independently of j) we can again solve this equation for d^A by contraction and obtain 
recursively estimates on the derivatives up to order k with the right dependence in 9. The 
first part of the Lemma is then proven on the time interval [0, T] . 

We now observe that at time T the estimate is in fact better. Indeed what we said 
before is also true for the time interval [0, T/2] . We can now write another integral equation, 
namely for T/2 < r < T 

A(t, • ) = e^~ T ^G {T _ T/2) * A(T/2, ■ ) 
+ / dsG T - s * (a(p + ij)A(s, -)-(l + ip)\A(s, -)\ 2 A(s, •)) , 

JT/2 

and for the gradient 

VA(r, • ) = e «(-T/2) VG(T _ T/2) * A ( T /2, ■ ) + f VG r _ s * (a(p + i^d^A^ ■ ) 

JT/2 

-(1+0)\A( 8 , -)\ 2 A(s, •)) da, 

We then deduce from the integrability in t = of the L 1 norm (in x) of VG t an estimate 
on the L°° norm of the gradient of A on the time interval ]T/2, T] which is independent of 
6. Similarly we can control recursively higher order derivatives by writing down integral 
equations starting at various increasing times between T/2 and T. 

In other words, even if we started at time zero with a number 9 which was quite 
large, after a time of order at most one we get an estimate where the function and all 
its derivatives up to order k are of order unity. The result for all times then follows 
immediately from previous results^ 2 ). Finally the result on the time derivatives follows at 
once from the result on the space derivatives using equation (A.l). 

We now need to prove that when a is small, if at the initial time A is of order unity, 
then it will decrease. Note that if a = (3 = this follows easily from the maximum 
principle applied to the square of the modulus of A. In the general case, we will use a 
method of local energy estimate^ 3 ' 2 ). 
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Lemma A. 2. Given the numbers a, 13, p, S, and K as before, there is a number 6 > 
and a number cr > such that if a G [0, <7 ] and t > 0, we have 

\\A(t, .)|U~ <S- x /2. 



Note that we could get a bound of order a/cF but after some larger time which depends 
on cr. One can also derive similar estimates for the derivatives. 

Proof. The proof is essentially based on a local energy estimate^ 2 ). First of all we can 
get rid of the term i^aA by considering the function e~ %iat A{t, x). We introduce then a 
cut-off function depending on a positive parameter e 



<p(x) 



e d 



(l + e 2 \x\ 2 ) d 

where d = 1 or 2 is the dimension. One then introduces the function of time 



= J cp(x)\A(t,x)\ 2 dx . 



It is easy to verify using integration by parts that since A is bounded, this function satisfies 
^vjr(t) < - 2 J cp(x)\VA(t,x)\ 2 + 2ap J <p(x)\A(t, x)\ 2 

+2(1 + |a|) j \Vy{x)\\VA{t,x)\\A(t,x\-2 J ip(x)\A(t,x)\ 4 . 

Completing the squares and using that |V0|/0 < 0(1 )e we get for some number a > 1 
independent of e 

±^(t)<+2(ap + ae 2 ) j <p(x)\A(t, x)\ 2 - 2 J \p(x)\A(t, x)\ 4 . 

Using the last term on the right hand side to control the first one, we finally get 

j^it) < -(ap + at 2 ) J <p(x)\A(t, x)\ 2 + 0(l)(ap + ae 2 ) 2 . 

This implies since ^(0) = 0(1) that beyond a time of order (ap + ae 2 ) -1 we have 

< 0(l)(ap + ae 2 ) . 

Let now t be larger than the number U of the previous Lemma. Then we have 

\\VA(t,-)\\ L ~<0(l). 

From the previous estimate we also know that the average of \A\ 2 on any square (segment) 
of size e _1 is bounded by 0(l)(ap + ae 2 ). Given S, this implies that in dimension one we 
can choose do and e small enough such that the estimate holds. In dimension 2, one needs 
a similar energy estimate for the gradient to apply an adequate Sobolev inequality. We 
refer to the literature for analogous estimates. 
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Appendix B: Proof of Lemmas II. 4 and II. 5. 

In this appendix we will always assume that the parameter r\ is positive and smaller than 
one. We will first give an integral representation of the kernels of the operators L l and 
L f Rs } g. These operators are convolution operators with a sequence and we will estimate 
the I 1 norm of this sequence. The operator L* is the convolution with the sequence V? 
given by 

K = (27T)- 1 / *e-^(, + IM)^ (Bl) 
Jo 

and the operator L l Rs^ is the convolution with the sequence M 1 

m\ = fj e-^(v + my (i - * {^f) - * i^f)) <%> m 

since 2n/59 is larger than twice the diameter of the support of ip. 

First of all, since / is C 2 , we can find a positive number a < min(u;, n — ui) such that 
on the intervals [±u — a, ±uj + a] we have \l"(<p) + D\ < D/2. Moreover, there is a positive 
constant B < 1 such that 

sup \l((p)\<B. 

ip$]±u} — a,±uj-\-a[ 

In the integral (Bl) we separate the domain of integration in the union of the two 
intervals [±u — a, ±uj + a] and their complements. Using integration by parts, one checks 
easily that the contribution of the complement to the integral is bounded in Z 1 norm 
by 0(1) (1 + t) 2 (i] + By. For the contribution of each interval we obtain by standard 
arguments^ 9 ) 

Itt)- 1 ( " " e-^irj + li^Ydip < 0(l)e tr >t-^ 2 . 

•J znuj — a 

On the other hand, integrating twice by parts we get for t > 2 and n ^ 

n^Vi = -(27T)" 1 e~ in * + l(<p)) + t(t - 1)1' {^f) ( V + l{<p))*-*&p 

and by similar estimates as before, for n ^ and t > 2 

IVH < C(l)n- 2 t 1 / 2 e 7 '* . 
Combining the two estimates we get 

IVH < C »(1)(1 + n 2 1 <t)-H- y /2 ! e"* . 

which implies immediately 

\\V*\\ii <0(l)e"*, 

and this proves Lemma II. 4. 

We now come to the proof of Lemma II. 5. First of all, if t < (56)~ 2 : the estimate 
follows at once from Lemma II. 4 and Lemma II. 1 (as long as 5 2 6 2 > r\ which is always 
assumed). From now on we will assume t > (59)~ 2 . 

It follows by estimates similar to the previous ones that 

\M f n \ < 0(l)(l + n7t)-Y/ 2 e-^ 2 * 2 / 2 -^ . 
from which Lemma II. 5 follows. 
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Appendix C. Estimation of oscillating sums. 

Lemma C.l. There is a constant Cg > 0, a constant a > and a constant 1 > rjo > 
suci that if 77 G [0, 770] , if 5, 0, B and A satisfy the hypothesis of Lemma II. 3, if the sequence 
/* is defined by 

f n = A 3 (S 2 t,nSy/2/D), 
then for any integers s and t we have 

\\L s e^f \\ioo < C 9 (59e vs + e~ as ) . 



Proof. From hypothesis H3 it follows that there is a number T > such that \4u 
(mod 27r) I > T and \2u> (mod 2tt)\ > T. We can now write for the kernel of the operator 

L s (n,m) = ±-( eW n ~ m \n + l{<j>)Y dcj> = ±- [ e*<»- m >fa + Z(0))'d0 



+ 



-/ 



D i4>{n—i 



|<T 



>{r] + l{(t>)) s d(t> + 



27r J|^+o.|<T 



|0±^|>T 



Using techniques analogous to those of the previous appendix, it is easy to show that the 
first integral is the kernel of an operator with norm 0(l)e~ as in Z°° for some constant 
a > which can be chosen independent of 77 if 770 > is small enough. The second and 
third integrals are estimated by similar methods of "summation by parts" and we will only 
treat the second one. 
We observe that 



-u>\<r 



YjU™ /m+i) 27r / 



^ n - m)+3lurn (r] + l((j))) s d(j) 



e iHn-m)+3iu>m^ v + 
]_ (ii(<t>— 3ui) 



'\<j}-Lo\<r 

From the properties of the function A (see Lemma A.l) we have 

\f* m - f* m+ i\ < o(i)se . 

Using again stationary phase methods, it follows as in the previous section that 



1 

27 



'\<j>-u\<T 



e —i(4> — 3uj) 



< 



0(l)s 



-l/2 pV s 



1 + (n — m) 2 / s 



and the Lemma follows. 
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Lemma C.2. Let (q n ,m) be such that 

^2(1 + \n\ + \m\) k+2 \q n , m \ < oo 



and q(4>i, 02) = on a neighborhood of the four points (±u;, ±u;). Then there is a number 
5 > such that if A is constructed as in III. 2 with an A(x) such that there is a constant 
M and a number k > such that 



sup 0'||dM|| < M , 

0<j<k 



then 



sup 

r 



^ ^ q-i — n,i — mA n A m 

n,m 



knk 



Proof. The proof again essentially mimics the proof of the stationary phase estimates, 
namely it relies on an integration by parts. Let F(x) = (3g)~ 1 / 2 A(5 2 t, ^/2/D(x + SVt)), 
we have to consider four terms which are essentially of the same form. We will only treat 
one of them in details, the others are estimated similarly. We have 

Y,Qr-n,r-me l ^ n+m)) F{5n)F{5m) = 

n,m 

e 2iu,r Y^F{5n)F{8m)-^ f ?(0i, fae^-^-^+fa-uKr-m))^ ^ _ 

n,m 

Let w be a positive number such that q(4>i, 02 ) = on the four sets [e\u — w, e\u: + w\ x 
[e 2 u; — w, e 2 w + w\ where e\ = ±1 and e<i = ±1 . We now decompose the integration 
domain into a finite union of domains, in each of which 0i or 02 (but may-be not both) is 
at a distance at least w of ±u. Let us call V such a domain, and assume for definiteness 
that on V we have 

inf |0i ± a; I > w . 

It is then easy to verify that 

J2 F (^)F(5m)-^ [ g(0i, 2 )e l(( ^- w)(r - n)+( ^- w)(r - m)) #i# 2 = 

Y^F{8{r-n))F{8{r-m))^ J qfa, 2 )e*«^-^+^-^ m ^0 1 # 2 = 

n,m ^> 

J2(F(8(r-n))-F(8(r-n + l))F(8m)^ f J^e^^-^^. 

n,m ^* 

From our hypothesis, we have 

\F(S(r - n)) - F(5(r - n + 1))| < C>(1)56> . 

The result follows by several applications of the argument, and controlling the convergent 
of the ensuing sum by the summability properties of (q n ,m)- 
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Lemma C.3. There is a number b > such that if < r\ < 1, if tie sequence (q n ,m) is 
such that 

+ |n| + |m|) 2 | 

and q(4>i, fa) = on a neighborhood of 4>i + 4>2 = then for any t > we have 

||£ t( ?(lZ«OI|z°° < C > (l)e — bt ||^;||ic« ||w;j| Z oo . 



Proof. We have easily the formula 

The result follows at once from the fact that on the support of q we have (pi + 4>2 7^ ±<^- 
The necessary summability results are obtained by integration by parts as in the previous 
Lemma. 
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